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Abstract 

A simple method to compute QED bound state properties is pre- 
sented, in which binding energy effects are treated non-perturbatively. 
It is shown that to take the effects of all ladder Coulomb photon ex- 
changes into account, one can simply perform the derivative of stan- 
dard QED amplitudes with respect to the external momentum. For 
example, the derivative of the light- by- light scattering amplitude gives 
an amplitude for orthopositronium decay to three photons where any 
number of Coulomb photon exchanges between the e+e^ is included. 

Various applications are presented. From them, it is shown that 
binding energy must be treated non-perturbatively in order to preserve 
the analyticity of positronium decay amplitudes. 

Interesting perspectives for quarkonium physics are briefly sketched. 
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1 Introduction 



The properties of positronium, the bound state made of an electron and 
a positron, provide some of the most precise tests of QED [0. Both the 
experimental and theoretical considerations have reach a very high level 
of precision, requiring for the latter the computations of many-loop diagrams 
references cited there). 

In the present note, we will address one particular aspect of the current 
QED bound state models, namely the factorization between the bound state 
dynamics and the decay processes |^] . As explained in |^ , the present study 
is motivated by the recurrent contradiction between standard positronium 
models and Low's theorem |^. Our goal is to show that binding energy 
effects cannot be treated perturbatively if sensible theoretical results are 
sought at the present level of precision. The central result of the present note 
is a simple method that allows such an exact non-perturbative treatment of 
binding energy effects. The lowest order results we shall obtain sum infinite 
class of diagrams, usually expanded perturbatively. Also, our method is 
formally relativistic and has a correct analytical behavior. 

In section 2, our simple method is introduced by a computation of the 
rate for parapositronium to two photons, and then demonstrated in general. 
In the next sections, we apply our result to various decay modes chosen to 
illustrate some specific aspects. 

First, to show explicitly the contradiction between Low's theorem and 
binding energy perturbative expansions, we consider the paradimuonium 
(p-Dm, a bound /^^/U^ [10|) decay to je~^e~ . Then, interesting connections 
between bound state decay amplitudes and the photon vacuum polarization 
are exemplified by the orthodimuonium decay o-Dm e^e~ . 

The fourth section is devoted to the interesting decay mode orthopositro- 
nium to three photons. We obtain the differential spectrum and decay rate 
for any value of the binding energy, from the amplitude for light-by-light 
scattering. We discover that both the spectrum and rate are highly sen- 
sitive to the binding energy when it is approaching zero. The implication 
of this result is the well-known slow convergence of the orthopositronium 
radiative correction perturbative series. 

In the fifth section, we present two extensions of our method. The appli- 
cation to orthopositronium nth radial excitations o-Ps (nS) is introduced by 
recomputing the decay rate — > jo-Ps (nS) Finally, it is shown how 
the positronium hyperfine splitting, i.e. the mass difference of the ortho- and 
parapositronium, can be obtained from vacuum polarization graphs through 
finite mass renormalizations. 
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2 Lowest Order Decay Amplitudes 



Positronium amplitudes are built as loop amplitudes: the positronium cou- 
ples to a virtual e^e~ loop, to which a given number of photons are attached. 
The coupling of the positronium to its constituents is essentially determined 
by the positronium wavefunction. In the following, we will take a form factor 
inspired from the Schrodinger wavefunction, i.e. a form factor containing 
the effects of the exchange of Coulomb photons among the constituents (in 
the ladder approximation) |12|. 



A very simple method can be used to compute such loop amplitudes: any 
lowest order loop amplitude with a Coulomb form factor is the derivative with 
respect to the positronium mass of the corresponding point-like amplitude. By 
point-like amplitude is meant the loop amplitude obtained by replacing the 
complicated Coulomb form factor for the bound state by a constant from 
factor. Symbolically, for parapositronium (orthopositronium) decay to an 
even (odd) number of real or virtual photons, the amplitude is 




oc 



d 

Qp2 




where (p represents the Schrodinger wavefunction, i.e. 





and dashed photon lines stand for Coulomb photons. 

In the next subsection, we give an example to illustrate explicitly how 
the method works, and in the following one, we present the general demon- 
stration. 

2.1 An Example: Parapositronium to two photons 

The point-like amplitude for the decay p-Ps 77 is obtained by replacing 
the Coulomb form factor by an elementary, point particle, with a pseu- 
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doscalar 75 coupling to the electron current (with unit coupling constant) 




plus the crossed process. The point bound state amplitude is fully relativis- 
tic, and standard loop integration techniques can be used. Explicitly, the 
point-like amplitude is (the subscript 75 is a reminder for point-like): 



M,, [p-Ps ^ 77) = ie^ I -^Tr e, (h) {I 

J (27r) 

T^" = 75 ^ 



1 



1 1 



+ 75- 



j/2 - m 
1 



i+]/2 -m 
Carrying the trace, we readily obtain 



m 



where the dimensionless loop integral form factor is 



J/i — m 



M2 



7F 



4m 



p2 



with F [a] =2 arctan'^ (a — 1) 



-1/2 



From the amplitude, the decay width is 

^2 2m2 vrMa^ f ij'^ 



{p-Ps 77) = IQira' 



m 



256 



2 M\ 
arctan — 



(2) 



27; 



where 7^ = — M'^ jA is related to the binding energy Eb = M — 2m. 

To get the physical positronium decay amplitude and rate, simply replace 
in (^) the loop form factor X^j. (Af^) by its derivative 



Icoui (M2) = (327rC<Ao7) ^^^5 {M^) 



(3) 
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where (po is the ^-wave fundamental state Schrodinger wavefunction at zero 
separation, and C = \fM jra is obtained by matching the static hmit (7 — > 0) 



with the weU-known lowest order result T:p.pg = ma 



72 0. This gives 



'-Coul 



M 



2 M 
— arctan — 
vr 27 



The factor in square brackets is equal to 1 in the limit 7 
|(/)oP = a^m^/SiT, the decay rate into two-photon is 



0. Using 



r {p-Ps — > 77) 





2 M 




— arctan — 


TT 27 



(4) 



The result @ contains the effects of Coulomb interactions among the con- 

M^/i the 



m 



stituents, at all orders in a. Indeed, using 7 
form factor can be expanded as (see the discussion in 



r {p-Ps — > 77) 



l-2- + 0(a' 

IT 



^1 - 0.6370 + 0(02)) 



In other words, the binding energy effects included in our lowest order com- 
putation already account for a great deal of the relativistic and radiative 
corrections as presented in the literature P]: 



p-Ps 



(1 - STp_Ps) 



(5) 



where 
6rp.ps 



Ap- + 2a^ In - + Bp— 

TT a TT^ 

A„ = 5- 7r74 ^ 2.5326 

with ^'^ r 1 A /nn\ 

Bp = 5.14(30) 



Numerically, we can write (|5|) as 

a^m Aw? 



3a^ , n 1 ^ , 1 ^ oP- 
\t? - + C„ — In - +5a 



2tt 



a 



Cp - 
64^ 



TT a 
-7.919(1) 
0.274(1) 



p-Ps 



2 M\ 
— arctan — 
TT 2j J 



2 

1 (1 



2 M2 

with the same series 0, but with the reduced coefficients 



^^p-Ps) 



A'p 
B'p 
Cp 



0.5326 

0.607 

-3.919 



TT^ 



(6) 
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This last form is very interesting because binding energy corrections (i.e. 
7-dependent) are singled out, while the remaining radiative corrections are 
much reduced. This means that one could, at least in principle, express the 
decay rate as non-perturbative binding energy corrections times a rapidly 
converging perturbation series of radiative corrections. This is exactly what 
we have achieved to order a. 



2.2 Generalization 

The proof of is straightforward using the language of dispersion relations. 
Indeed, the point-like loop amplitude can be computed from its imaginary 
part 

ImT^, (P2) ^ i,^M^, {p-Ps ^ 
using an unsubstracted dispersion relation with s = P^: 

r,, (M^) = Re {M') = - J^^^ Im T,, (s) (7) 

(T (M^) = ReT(M^) because M < 2m). The Schrodinger form factor 
accounting for the non-trivial coupling of the bound state to its constituent 
is of the form 



FB(0,q) = C(Ao.F(q2) (q2+72) with (q 



2^ 



87r7 



(q2 + 72 



,2 



where one can recognize (j) (q^) = 4>oJ' (q^) as the fundamental (n = 1) S- 
wave Schrodinger momentum space wavefunction. When expressed in terms 
of the dispersion relation variable, this form factor is only a function of s, 
the initial energy: 

Fs is) = C^o- ''"^ 



M2 

The core of the derivative approach emerges from the observation that in- 
serting Fb in is equivalent to taking the derivative with respect to 

1 /■+°° ds 

1 r+°° ds 

= (32vrC<Ao7) - / 77^Im^75 

^ J4m2 [S - M^Y 

= (327rC<Ao7) ^r,, (m2) (8) 
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which is the desired result. The case of other parapositronium decay chan- 
nels, or orthopositronium decay modes is similarly treated (simply replace 
7^ by ^ with the orthopositronium polarization vector). 

As a first remark, we repeat here the conclusion of Q, which is that 
standard positronium result are recovered from the above dispersion integral 
(0) provided only the vertical cuts are taken into account in the imaginary 
part. This proves that some contributions are missed in those computations. 
For instance, in the case of p-Ps — > 77, there is only the vertical cut, but 
this is not the case in general: if one of the photon is virtual, or if there is 
three or more photons in the final state, oblique cuts contribute (see next 
sections). 

The second remark concerns the insertion of the form factor, not in 
Feynman amplitudes, but directly into the dispersion integral. These two 
approaches are equivalent only in the punctual case (consider the momentum 
flow through the form factor in each case). Working at the level of dispersion 
relations is much more in the spirit of the Bethe-Salpeter equation. Indeed, 
the wavefunction is extracted from the four-point Green's function, i.e. in 
configurations with an off-shell bound state (above threshold), and on-shell 
constituents: 




The picture shows that both Bethe-Salpeter and the dispersion integral (0) 
make use of the form factor Fb with the same kinematical configuration 
(diagrams with bremsstrahlung radiation off the electron lines are treated 
similarly) . 

3 Application to decay rate computations 

We will now review some applications of the result (|8|) . We leave the inter- 
esting case of o-Ps 777 to the next section. 

3.1 Paradimuonium and Low's Theorem 

The paradimuonium decay p-Dm — > 7e'^e~ is the simplest QED bound 
state decay process where Low's theorem implications can be illustrated [^, 
[§, [^. Following the same steps as for p-Ps — > 77, we first consider the 
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point-like amplitude 




(plus the crossed diagram) with the result 

1 — ae 



■ 75 



(p-Dm 



e e 



-7) 



3 W 



with p{x^,ae 



1 — X. 



J dx.y |Xyg {M"^, x^) p {x^, tte) 


[2 (1 - X^) + Oe 



(9) 



-i 1 








"4m2 1 






-F 




(47r)^ 


_M2 


MM - 



where m is the muon mass, M the dimuonium the reduced photon 

energy 2E^/M, = 4mg/M-^, me the electron mass, and with the function 
F defined in (|l]). In the limit x^ 0, the spectrum dV^^/dx^ goes to zero 
as x^ as predicted by Low's theorem (the amplitude behaves and an 

additional factor x-y comes from phase-space). 

Taking the derivative of Xy^ to get the corresponding Coulomb form 
factor, we find 



ICoul (M^,x, 



(327rC(/.o7) 
.Ccko 1 



where 



M x- 
4m 



2 M 

— arctan arctan 

vr 27 vrM 



1 



-1/2 



,M2(1 -x^) 

The decay rate is then obtained by replacing I.y^ by Icoui in 



r (p-Dm — > e~^e 7) 



a^m 4m2 
67r M2 



l-a. 



dx^ 



2 M 47^^, 

— arctan arctan 

vr 27 vrM 
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It is very instructive to analyze in some details this result. Consider the 
limit for the form factor: 

Wmlcoui {M\x^) = -i^ (10) 

In that limit, the standard result for the decay rate is recovered 



i [p-Um — > e e 7) = / ax. 

67r 7o 



7 „2 



However, the differential rate dF/dx-yhas a wrong behavior when — > 0. 
The spectrum is linear (in x-y) in the limit 7 — > 0, in contradiction with 
Low's theorem. Therefore, it appears that, contrary to the two real photon 
case, the limit 7 ^ is far from smooth. It is inconsistent to consider both 
the soft-photon limit and the on-shell limit simultaneously. Explicitly, the 
incompatibility of the two limits is obvious if ^ is taken first 

h^lcoul (M , X,) = + arctan - J (11) 

= /M2 1 _ 

* M V872 2 SMvr 

Mathematically, what these considerations show is that the limit 7^0 
does not exist at x^ = 0. Consequences of this result were discussed in [^, 
but it should be clear that one of the main virtue of our approach is its 
non-perturbative treatment of 7, leading to correct photon spectra. 

To close this section, we just mention that the second term oiXcoui (-^^j ^7) 
can be traced to the presence of oblique cuts in the imaginary part (i.e. pro- 
cesses like p-Dm fi'^fi^^ times fj-^fi~ (7) — > e^e~ (7), see Q, Q). Those 
are the processes neglected in standard approaches, which is an inconsistent 
approximation since it is the sum of the two terms of Icoui {M'^,Xy) that 
enforce Low's theorem. This can be seen by plotting 

J(x„7) = l- ^"^'^"^:7^ 
arctan 
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for M = 1 and various values of 7 




0.2 0.4 0.6 0.8 1 



As the picture shows, even if 7/M << 1 for QED bound states, the hmit 
7 ^ is not to be taken because it is singular. 

3.2 Orthodimuonium and Photon Vacuum Polarization 

We now consider the decay o-Dm — > 7* — > e~^e~ 



o-Dm 




The point-like amplitude is easily obtained in terms of the (divergent) pho- 
ton vacuum polarization function 



M^. {o-Dm ^ e+e-) = (P) {P^g^'' - P^P'') H^m (P^) i_ [u (p) [p') ] 
where 



n^M (p2) 



47r2 



(12) 



with Q = P"^ /m?, m the muon mass and D = 2/e — "fEuier + log Air/i'^/Tv? in 
dimensional regularization. 
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The Coulomb form factor will be obtained from the derivative of the 
vacuum polarization, with respect to = M"^: 



^Coui (m2) = (327rC70o7) 



d 



eC4>o 

M 
eC4)o 

M 



7- 



3M3 



SvrM 



'M2 



r arctan — 

27 y 



128 



(13) 



This has the effect of removing the divergence, as it should. The decay rate 
is 



r (^o-Dm e^e ) 



a M2 + 2ml 



M 



VI - ae In 



Coul 



~6~ 



SvrM 



'M2 



To leading order, we recover the standard result T (o-Dm e^e~) ~ 

when Oe << 1. For electromagnetic bound states, the binding energy is 

related to the fine structure constant, hence the corrections can be cast into 



r (o-Dm 



e e 



6 



1 + 



1 + 



16 a / 2 64 

3 TT 



+ .. 



1 - 1.70a + 1.720^ + ...) 



Compared to the binding energy corrections to the parapositronium two- 
photon decay rate, the present corrections are much bigger. Note also that 
the binding energy correction obtained here, at order a, again account for 
a great deal of the correction obtained using the standard approach (we 
consider only part of the total O (a) correction, see Eq. 47 in [^), which 
are 



^LO+NLO / 



e e 



a^m 



1 



4^ 



4 Orthopositronium Decay 



We now apply our method to the orthopositronium decay to three photons. 
This is a very interesting decay process. As discussed in Q, the lowest order 
basis chosen in standard computations, namely the Ore-Powell amplitude 
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|13|, is in contradiction with Low's theorem. What we will now show is that 
taking binding energy into account, i.e. integrating the effects of all the 
Coulomb photon exchanges at lowest order through the appropriate form 
factor, is necessary in order to get a correct spectrum. Also, we will again 
find that the bulk of the radiative corrections is accounted for already at 
our lowest order. 



4.1 Point-like Amplitude 

The point-like amplitude is the standard light-by-light box diagram 



o-Ps 




plus five other ordering of the photon insertions. The amplitude can be 
found in many places The tensor G'^^'^'^'^^ describing the transition 

from an off-shell photon to three on-shell photons of helicity states A1A2A3 
is such that 

fQXi\2\3Q*\i\2X3,a\ ^ _^ (123) + R (213) + R (312)] (14) 

where i?(123) = R{xi,X2,x^,a) {xi is the reduced energy of the photon i 
and a = Am? /M'^). The R (ijk) are given in terms of individual dimension- 
less helicity amplitudes as 



R{123) 



+ 



i|i?it+(123) 

xi 

X2X3 (1 - Xl) 
(1-X2)(1-X3) 



+ 



(2) 



(123) 



E_ 



(1) 



(213) 



+ — 

xi 



E 



(1) 



Xf (1 - Xl) 



1 



X2 



-E^A+ (123) 



+++(123) + i?+ 
1 ^(1) 



(1) 



(132) 
2 



(15) 
2 



1 - X3 



-E\ 



(132) 



The helicity amplitudes -£^±4.+ are complicated functions of Xi and a, and 
we do not reproduce them here (p, IQ). The decay rate of a point-like 
vector positronium to three photons is then 

1 



(o-Ps 777) 



--^M / dxidx2 [R (123) + R (213) + R (312)] 
9 z°7r^ J 
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with the three-body phase space written in terms of reduced photon energies 
as / d<^3 = AP / dxidx2/2'^7r^ §. 

Of special interest is the behavior of the low-energy end of the differential 
rate. For various values of the ratio a = 4m^/M^ (i.e. of the binding energy 
Eb = M — 2m), the photon spectrum, normalized to the total rate is 



1 dE 



100 




Sufficiently close to zero, the behavior is always like required by 

Low's theorem |^]. As a increases, the behavior is getting more and more 
pronounced. In the limit a ^ oo, the normalized spectrum is 



1 fir. 



a— >oo 



r^M dxi 



17 



^3 

Xi 



343 
~3~ 



973 

207x1 + —X 



which is, as expected, the spectrum obtained from the Euler-Heisenberg 
effective theory [IS] 



^E-H 



a 



90m4 



4.2 Coulomb Form factor and Ore-Powell spectrum 

To insert the Schrodinger wavefunction form factor, we define modified he- 
licity amplitudes according to (S) 



(327rC0o7) 



9 pin) 



{ijk,M^) 
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The resulting decay amplitude, and decay rate are constructed as in the 
punctual case, and we reach 



r {o-Ps 777) 

with TZ = 



4m ^ 



7^ 
2^ 



dxidx2 {Rcoul (123) + Rcoui (213) + Rcoui (312)) 

Rcoui is given by (jl^) with E^^j^ ^^^^^ in place of • Integrating over 

X2, we get the photon spectrum for various values of a 



rdx 

2 

1.75 
1.5 

1.25 
1 

0.75 
0.5 
0.25 




(Ore-PoweU) 

a-^1 

a=l + 10^ 

a=1.01 

a=1.05 
a=l.l 
a=1.5 

a— >oo 
(Euler-Heisenberg) 



0.2 



0.4 



0.6 



0.8 



1 



X, 



(16) 



As long as a 7^ 1, the spectrum behavior is in close to zero, i.e. roughly 
in the range xi G [0, a — 1] . The properties of this spectrum as a varies are 
completely similar to that p-Dm je~^e~, and one can show that the 
limits 7 — > and xi — > are again incompatible. 

In the context of dispersion relations, the contributions originating in 
the oblique cuts (processes like o-Ps e^e~^ times e+e" (7) 77 (7) in 
the imaginary part) are essential to maintain a physical spectrum (i.e. in 
agreement with Low's theorem). Again, those oblique cuts are neglected in 
standard approaches. Here, they are automatically accounted for since they 
are included in the point-like result. 

Note also that the spectrum for a point-like bound state at threshold 
corresponds to the spectrum for a ~ 1.05 in the Coulomb form factor case. 
The bound state decay spectra for a < 1.05 are unattainable in the point 
bound state case. 
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Concerning the total integrated rate, we recover the Ore-Poweh result 
at threshold (a = 1) 



r {o-Ps 777) = a^m- 



17.16 



4.75 X 10"^^ MeV = a'^m— (vr^ - 9) 

yvr ^ ' 



'<dTT V 27r2 ^ 

As a increases, the total rate quickly decreases (remember a = \r^'^ jM'^ + 1 

11 

lie 

ie-9 
0.8 

0.6 

0.4 

0.2 



.001 .002 .003 .004 .005 



By a numerical analysis, we found 



= (tt^ - 9) (l - 15.412^ + 122-^ - 889-^ + 1.92 10^ ^ - ..\ 

For orthopositronium, we can express 7 in terms of a, and we find the 
binding energy corrections to the total rate 



r (o-Ps — > 777) = a^m- 



2 7r'-9 



9tt 



1 

a a a 

1 - 12.1- + 80.2^ - 502^ + ... 



(17) 



This series is to be compared to the one presented in the literature Q, which 
is 

2(7r2 



To-ps = a m- 



97r 



1 



Ao ^ In - + - — In^ - 



CK 1 a 
-Co — In - + ^5^- 

TT a TT 



with coefficients 



Ao = 10.286606 (10) Co = 5.517 (1) 
5o = 44.52 (26) (55^ = 0.19(1) 

Again, the exchanges of Coulomb photon appear responsible for the bulk of 
the radiative corrections at order a. Note that the origin of the slowness 
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in the convergence of the corrections to Tq-Ps is clearly identified as coming 
from the perturbative expansion of binding energy effects (see (|l7|)). In 
other words, the four-point fermionic loop with a Coulomb form factor is 
not well-behaved for 7 = 0, and that limit appears as an inappropriate basis 
for perturbation theory. 



5 Application to other processes 

In this final section, we present two extensions of the method. First, we 
accommodate the derivative formula (P) for radial excitations, by consider- 
ing the pion decay vr'^ — > jo-Ps. For the second, we compute the hyperfine 
splitting using the language of mass renormalization. 



5.1 Pion Decay to Orthopositronium 
and Radial Excitations 

The technique of taking the derivative of point-like amplitude can be ex- 
tended to other spherically symmetric wavefunctions. The decay ^ j o- 
Ps is a good example [0] 




To get a sensible theoretical prediction, one must sum the decay rates over 
the infinite tower of radial excitations o-Ps {nS) in the final state {o-Ps {nS) 
means nth radial excitation of the 5- wave J = 1 positronium state). From 
the standard Schrodinger wavefunction for hydrogen states, we can write a 
general expression for the form factor for nS radial excitations as 



-CouLn 



[Ml) = {^2tiCc^^ 



iFi 



(19) 



with 



ma 
Svrn^ 



7n = - Ml/ A 
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Where we have denoted the mass of the o-Ps (nS) state. The hyperge- 
ometric functions are essentially the well-known Laguerre polynomials ((|T9|) 
is equally valid for parapositronium). 

For the case at hand, is the photon vacuum polarization (|l2|), where 
m refers now to the electron mass. The pion decay rate into orthopositron- 
ium states can be written as 

T (tt^ ^ o-Ps (nS)) / 9m2 / M^\^ \ 

RoPsinS)= r(vr0^77) -''^"\^Coul,n{M'n)\ " 4j + ^ 

The M^/rnji^^ ~ 10^^ corrections arise from the form factor for pion to 
two photons. The mass ratio M^/m^ ~ 10~^ is also negligible compared 
to binding energy corrections, to which we now turn. Up to corrections of 
order 7^, we can write using 7„ « ma/2n 

(the An are the numerical coefficients found by expanding Ilcoui,n (^n))- 
Summing over n, we find 

4 

E ^o?sinS) = yC (3) (1 - 1.66a) « 1.684 x IQ-^ 

n 

where ~ (1.418 x 10^^) is obtained from the contribution of the o- 



Ps (IS) only. For comparison, [11| found the radiative correction to be 

3d branching frac 

(1.9 ±0.3) X 10" 



(1 — 0.92a). The experimentally quoted branching fraction is [16| 
r (tt ^ -fPs) 



r(7r0^77) 

so the agreement is good. 

What this little exercise shows is the power of our method as a mean to 
partially compute higher order corrections. More importantly, it is by now 
apparent that other, more complicated wavefunctions can easily be accom- 
modated for. Any wavefunction that can be expressed, even approximately, 
from derivatives of the Coulomb one can fit in our scheme. This may open 
the way to many applications in QCD, as we will comment in the conclusion. 

5.2 Hyperfine Splitting 

and Mass Renormalization 

We implement here the renormalization of the positronium state. The idea 
is to consider the bare positronium mass M as equal to twice the electron 
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mass m, and then to carry a (finite) mass renormalization. The diagrams 
contributing, at lowest order, to this mass shift will be obtained as the 
second derivatives of the vacuum polarization loops: 




Let us first obtain the mass renormalization equations, and then discuss the 
results numerically. 

5.2.1 Mass Renormalization 

For a pseudoscalar parapositronium state, the resummation of the Dyson 
series is trivial. The parapositronium propagator is then 



Of interest to us is the mass shift, defined from the pole of G {cp'^ 



Mlpara " = ^para {MI) ^ Mr,^,™ - = ^ L (20) 

For the orthopositronium, the transverse part of the bare propagator is 



G7 {i) 



- + ie 

The self-energy of the vector positronium will be of the form 
Proceeding with the resummation of the Dyson series, we end up with 



G^^ (g2) 



g2 (1 - Hortho (q^)) - M2 

Again, the propagator pole is at g2 = M^, hence 

Mlortho - Ml = Mlliortho {Ml) ^ Mn^ortho " = ^^^^^liortho {Mlortho) 

(21) 
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The two equations (20) and ( [2l| ) are self-consistent equations. By iden- 
tifying Mo = 2m and Mji^para{ortho) = Mpara{ortho) , they are of the form 

n Cm^ 

4m ^ 

Our approach wih be to generate corrections to )-,y plugging the 

non-relativistic result Eb = — ma^/4 into fP"-i~a(ortho) ^ 

5.2.2 Second Derivatives and Hyperfine Splitting 

The double derivatives of the vacuum polarization loops are 

1 1 / (9 \ ^ 

^paraiortho) {M^) = ^ (327rC(/.o7)' - ( ^^Hy^^) (M^) 

where a factor of 1/2 corrects for the factor 2 appearing in the derivative of 
1/ (s — M^) , while the other accounts for the double counting of Coulomb 
photon exchanges. 

The computations are straightforward. For parapositronium, the point- 
like quadratic divergence disappears and we find 

= M,^ra - 2m = "^°-]^'"-) 



1^ 

Mm 



472^^^arctanM ^ _ i' 



M=Mr, 



ma 



^ 9ma^ ma^ 115ma^ 



4 128 Gtt 8192 



+ ... 



Where we have used |0o| = m^a^/87r,7^ = m^ — Mp„^„/4 and Mpara = 
2m — ma'^/A in the right hand side. Proceeding similarly with the or- 
thopositronium renormalization, we consider the second derivative of the 
photon vacuum polarization function: 



E^nff ^ Mortho - 2m = ^^Uortho {Ml,,o) (22) 

|2 



mi 

M2 



.2 \ /on 2 \ arctan ^ o+i^^+ieo^ 



M=Mortho 



ma 27ma'^ 2ma^ ISOlma^ 
4 128 3^ 8192 "■ 
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At the order a^, there is also the annihilation diagram 




whose contribution is U'^^^t^o^'" [M^) = -i {P^g''^ - PT"") {Ucoui {P^)) 



Ann,fJ,U /- T, • / 7-.9 iiij n// /tt / t-,9.\\2 

ortho , . 

Using Hcoui {P^) given by (13), the dominant contribution in the limit 
7 — > is 

4 



^'■ortho 



E 



ortho, Ann 
B{2) 



Mprtho JT Ann 
o ortho 



{M ortho) 



ma 



+ .. 



Hence the lowest loop contributions to the hyperfine splitting are, to order 

,4. 



a 



AE. 



hf 



E 



•ortho 
B(l) 



E 



ortho, Ann 
B(2) 



^para 



ma 



17 

32 ^ 



ma'^ [0.5313 + 



This is to be compared to the result (see 
there) 



17 1 and references cited 



hf 



ma 



8 In 2 

77 + 



+ 



a 



a 



^■3 



ma 



7 a 

12 ~ ^ V9 ' 2 

5 2 1367 

-— TT InaH : — 

24 648 

--vr^ln^a + ( — ln2 
8 V 3 

[0.5833 - 0.3933a 



5197 9 / 
3456 V 144 



221 , 1 



ir'^ In a 



217 
"90" 
0.2083a^ In a 



2. 

+ 0{a'] 
0.3928a^ - 



In 2 



53 
32 



C(3) 



where the correction of order O (a'' In a) has been obtained very recently 

It is not surprising that our method does not reproduce exactly the above 
result, because we have neglected many diagrams (like electron self-energy 
insertions for example), and because we used the lowest order non-relativistic 
binding energy Eb = —ma^ /A and wavefunction |0o|^ = m^a^/87r as a basis. 
Taken individually, the corrections of O (a^) to the para- and orthopositro- 
nium masses are off by more than 50%. On the other hand, the difference 
between both corrections, giving the hyperfine splitting, is surprisingly good. 
Again, it seems that all the effects contained in the lowest order loop (i.e. 
all the ladder Coulomb photon exchanges) suffice to account for most of the 
radiative corrections. In conclusion, further studies of the application of our 
method to hyperfine splitting appear as necessary. 
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6 Conclusion 



We can characterize our method as a change of basis for perturbation the- 
ory: some binding effects are treated at aU orders in 7, at each order in 
a. With such non-perturbative results, the behavior of the usual expansion 
in the binding energy can be analyzed. We find that such expansions are 
problematic when soft-photons are present (i.e. when a photon's energy is 
of the order of 7). In other words, expansions in the binding energy are in- 
compatible with analyticity. In all cases, unphysical contributions violating 
analyticity are subleading, being of the order of 2m — M ~ How- 
ever small, such effects become important in view of the precision achieved 
in both experimental and theoretical descriptions. In addition, standard 
positronium models neglect some contributions (corresponding, in the con- 
text of dispersion relations, to oblique cuts contributions to the imaginary 
part of the loop diagrams). Therefore, they cannot be correct at order 
for decay processes [@]. 

In our view, obtaining QED bound state amplitudes from relativistic 
point-like amplitudes has many advantages. This may lead to simplifications 
in the computation of higher order corrections, especially concerning issues 
of gauge invariance and infrared divergence. For instance, two-loop graphs 
will introduce some O (q^) and O (q^ Ina) corrections (see for example |18]). 
We should point-out, however, that the technique of taking the derivative 
will not work as it stands for higher order diagrams, because of binding 
graphs, and some refinements will be necessary. 

As a final comment, the present method can lead to interesting advances 
in quarkonium physics. Of course, some knowledge of the quarkonium wave- 
functions will be necessary. No matter the form of these wavef unctions, it 
should be possible to expand them in the basis made of the Coulomb wave- 
function and its derivatives, to which the present method apply. This is 
left to further studies. For now, as a first result, it appear quite obvious 
that the Ore-Powell spectrum used as a lowest order approximation for the 
photon spectrum in the quarkonium inclusive decays into hadrons + photon 
is incorrect. If quarkonia were Coulomb bound state, the spectrum would 
be the curve of (|T6|) with a 1.2 ^ 1.5 (since Am^/My^j^-^ ^ 1.25 and 
Anrijj/Mj^^ 1-46), i.e. a spectrum which is suppressed at high energy 
compared to the Ore-Powell one. QCD effects should not change much our 
conclusion. We think that the suppression at high energy observed by CLEO 
]l9[ is, at least in part, a manifestation of binding energy effects, giving us 
confidence in our results. 
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